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1 Introduction 

In this paper we study the semiclassical limit of the Schrodinger equation. Under mild regularity 
assumptions on the potential U which include Born-Oppenheimer potential energy surfaces in 
molecular dynamics, we establish asymptotic validity of classical dynamics globally in space and 
time for "almost all" initial data, with respect to an appropriate reference measure on the space 
of initial data. In order to achieve this goal we study the flow in the space of measures induced 
by the continuity equation: we prove existence, uniqueness and stability properties of the flow 
in this infinite-dimensional space, in the same spirit of the theory developed in the case when 
the state space is Euclidean, starting from the seminal paper [13] (see also [1] and the Lecture 
Notes [2], |3]). 

As we said, we are concerned with the derivation of classical mechanics from quantum me- 
chanics, corresponding to the study of the asymptotic behaviour of solutions ip^(t,x) = V't (^) to 
the Schrodinger equation 

' iedtr^ = -4 A# + [/# = H,i;f, 

(1) 

as e — 7- 0. This problem has a long history (see e.g. |26]) and has been considered from a transport 
equation point of view in |22) and |19) and more recently in [7], in the context of molecular 
dynamics. In that context the standing assumptions on the initial conditions -00, e G H'^{W^;C) 
are: 

I |V'o,.|'dx = l, (2) 
sup / iH^ipQ^^l'^ dx < oo. (3) 

e>0 JR" 

The potential [7 in ([1]) is assumed to satisfy the standard Kato conditions U = Ub + Us with 

Us{x)= Yl y^xa-xp), VapeL\R^) + L'^iR^) (4) 

l<a<l3<M 

and 

Ub G L°°(R"), (5) 
VC/b G L°°(]R";R"). (6) 



Here n = 3M, x = [xi, . . . ,xm) G (M^)*'^ represent the positions of atomic nuclei. Under 
assumptions (|4]), ([5]) the operator Hgr is selfadjoint on L^(]R";C) with domain ff^(]R";C) and 
generates a unitary group in L^(]R"; C); hence Jjg„ \ipf\'^dx = 1 for all t G M, 1 1— s- tpf is continuous 
with values in if^(]R"'; C) and continuously differentiable with values in L^(]R"'; C). Prototypically, 
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U is the Born-Oppenheimer ground state potential energy surface of the molecule, that is to say 
Vai3{xa - xp) = ZaZp\xa - xp\~^, Z G N, \J\,{x) = inf spec iJ^f (x), where 



TV 



M 



Hel{x) = ^{-]^j:^r,-^Zo\ri- Xa\ ^) + ^ \ 



i=l 



a=l 



is the electronic Hamiltonian acting on the antisymmetric subspace of 



X Z2)^;C) and 



the Tj G M'^ are electronic position coordinates. For neutral or positively charged molecules 
(-^ < Sq=i Za), Zhislin's theorem (see |18) for a short proof) states that for all x, Ub{x) is an 
isolated eigenvalue of finite multiplicity of H(,£{x). 

In the study of this semiclassical limit difficulties arise on the one hand from the fact that 
Vf/ is unbounded (because of Coulomb singularities) and on the other hand from the fact that 
VC/ might be discontinuous even out of Coulomb singularities (because of possible eigenvalue 
crossings of the electronic Hamiltonian H^i). Fortunately, it turns out that these two difficulties 
can be dealt with separately. 

If we denote by b : M^" — )• M^" the autonomous divergence-free vector field b(x,p) : = 
(p, — V[/(x)), the Liouville equation describing classical dynamics is 



dtfJ-t + P ■ ^xfJ't - VC/ (x) ■ Vpfit = 0. 
If we denote by Ws : L^(M";C) — )■ L°°(]R" x M") the Wigner transform, namely 



W£'i/;{x,p) := 



1 



(27r)" 



Hx + ^y)fp{x - ^y)e 'P^dy, 



(7) 



(8) 



a calculation going back to Wigner himself (see for instance |22] or [7] for a detailed derivation) 
shows that Weipf solves in the sense of distributions the equation 



dtWei^t + p ■ V^Wei^t = #), 
where (§'ir{U,ip){x,p) is given by 



(9) 



4(C/,^)(x,p) := 



C/(x + |y)-C/(x-|y) 



^Pix + '-y)iP{x - '-y)e-'Pydy. (10) 



Adding and subtracting VC/(x) • y in the term in square brackets and using ye = iVpC 

an integration by parts gives ^'^{UjiJj) = VU{x) -X/pWeip + S'^{U,'ip), where S'^{U,'ip){x,p) is given 

by 



<(C/,V)(x,p) := 



{27ry 



U{x + ^y)-Uix-ly) 



-{VU{x),y) 



Hence, W^V't solves ([7| with an error term: 



i^ix+^y)i^ix - |y)e 'P^dy. 

(11) 

(12) 
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Heuristically, since the term in square brackets in (jlip tends to when U is differentiable, this 
suggests that the hmit of Wsipf should satisfy ([7|, and a first rigorous proof of this fact was given 
in |22) and |19) (see also |20)): basically, ignoring other global conditions on U, these results state 
that: 

(a) regularity of U ensures that limit points of W^i/jf as e J, exist and satisfy ([7|; 

(b) regularity of U ensures uniqueness of the limit, i.e. full convergence as e — t- 0. 

In (a), convergence of the Wigner transforms is understood in a natural dual space A' (see (|43p 
for the definition of A). 

In [7] we were able to achieve the existence of limit points even when Coulomb singularities 
and crossings are present, namely assuming only that Ub satisfies (O, (l6|, and (a) when Coulomb 
singularities but no crossings are present, namely assuming that Uf, G C^. If one wishes to 
improve (a) and (b), trying to prove a full convergence result as e ^ under weaker regularity 
assumptions on b (say S/U S W^'^ or S/U G BV out of Coulomb singularities), one faces 
the difficulty that the continuity equation is well posed only in good functional spaces like 
Lf{[0,T];L^ n L°^(M"')) (see [13], [1], [H]). On the other hand, in the study of semiclassical 
limits it is natural to consider families of wavefunctions ■0o,e in O whose Wigner transforms do 
concentrate as e ^ 0, for instance the semiclassical wave packets 



-na/2 



< a < 1 



(13) 



which satisfy lim^ W^V^Ce = ll'^o|li2'5(xo,po)- Here the limiting case a = 1 corresponds to con- 
centration in position only, lim^ WgV'Ce = ^xo ^ (2'?r)~"|J^0oP(" — Po)-^"; ^-nd the case a = 



9ol 



X (5pp. Here and below, 



yields concentration in momentum only, limg We^o,( 
{J~4>o)ip) = /jjn e~*^'^(/)o(x) (ix denotes the (standard, not scaled) Fourier transform. But, even 
in these cases there is a considerable difficulty in the analysis of (|10p , since the difference quotients 
of U have a limit only at ^"-a.e. point. 

For these initial conditions there is presumably no hope to achieve full convergence as e — )• 
for all (xo,Po)) since the limit problem is not well posed. However, in the same spirit of the 
theory of flows that we shall illustrate in the second part of the introduction, one may look at the 
family of solutions, indexed in the case of the initial conditions (|13p by {xq,pq), as a whole. More 
generally, we are considering a family of solutions tjjfyj to ([1]) indexed by a "random" parameter 
w £ W running in a probability space (W, J-, P), and achieve convergence "with probability one", 
using the theory developed in the first part of the paper, under the no-concentration in mean 
assumptions 



sup sup 

e>0 ieIR 



w 



w 



< oo, 



L°°( 



2n^ 



sup sup 

e>0 teK 



w 



L°°{R") 



< C{X) < oo VA > 0. 



(14) 



(15) 



Here Gg^"^ is the Gaussian kernel in M"^" with variance e/2. Under these assumptions and those 
on U given in Section 7.2, our full convergence result reads as follows: 



lim 

£4,0 



sup dj^>{Weiijt^,fj,{t,fi^))dF{w) = o vr>o 

Wt&[-T,T] 



(16) 
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(here d^/ is any bounded distance inducing the weak* topology in the unit ball of A') where 
IJb{t,iiw) is the flow in the space of probability measures at time t starting from and = 
lim^ WsipQ ^ depends only on the initial conditions. For instance, in the case of the initial 
conditions (fT3|) with \\(j)o\\2 = 1, indexed hy w = {xo,po), fi^ = 6^ and fi,{t,fiw) = 5x{t,w)i 
where X(t, w) is the flow in M^"" induced by {p, — Vf7). So, we may say that the flow of Wigner 
measures, thought of as elements of A', induced by the Schrodinger equation converges as e — )• 
to the flow in ^(M?^^ C A' induced by the Liouville equation, provided the initial conditions 
ensure (|14p and (jlSp . 

Of course one can question about the conditions (|14p and (|15p ; we show that both are implied 
by the uniform operator inequality (here p"^ is the orthogonal projection on ip) 

p^lw c]F{w) < CId with C independent of t, e. 

In turn, this latter property is propagated in time (i.e. if the inequality holds at t = it 
holds for all times), and it has a natural quantum mechanical interpretation. In addition, the 
uniform operator inequality is fulfilled by the classical family of initial data (|13p when P is a 
bounded probability density on M^'". These results indicate also that the no-concentration in 
mean conditions are not only technically convenient, but somehow natural. 

An alternative approach to the flow viewpoint advocated here for validating classical dynamics 
([7]) from quantum dynamics ([!]) would be to work with deterministic initial data, but restrict 
them to those giving rise to suitable bounds, in mean, on the projection operators p^"-^ . The 
problem of finding sufficient conditions to ensure these uniform bounds is studied in [T7]. Another 
related research direction is a finer analysis of the behaviour of solutions, in the spirit of |14) . |15) . 
However, this analyis is presently possible only for very particular cases of eigenvalue crossings. 

It is likely that our results can be applied to many more families of initial conditions, but this 
is not the goal of this paper. The proof of (|16p relies on several apriori and fine estimates and 
on the theoretical tools described in the second part of the introduction and announced in |Q. 
In particular we apply the stability properties of the I'-RLF in ^(M^"), see Theorem 15. 2[ to the 

Husimi transforms of ipf^, namely Wetpf^ * G^^\ Indeed, w*- convergence in A' of the Wigner 
transforms is equivalent, under extra tightness assumptions, to weak convergence in ^(R^") of 
the Husimi transforms. 

We leave aside further extensions analogous to those considered in |22], namely the conver- 
gence of density matrices p^, whose dynamics is described by iedtp^ = [H^, p^], and the nonlinear 
case when U = Uq * p being the position density of ^p (i.e. IV'P)- In connection with the first 
extension, notice that the action of the Wigner/Husimi transforms becomes linear, when seen at 
this level. 

Let us now describe the "flow" viewpoint first in finite-dimensional spaces, where by now the 
theory is well understood. Denoting by : M*^ — t- M*^, t S [0,T], the possibly time-dependent 
velocity field, the first basic idea is not to look for pointwise uniqueness statements, but rather 
to the family of solutions to the ODE as a whole. This leads to the concept of flow map X{t, x) 
associated to b, i.e. a map satisfying X{0,x) = x and X{t,x) = 7(t), where 7(0) = x and 

j{t) = bt{-f{t)) for ^i-a.e. t G (0,r). (17) 
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for .Sf^-a.e. x G M"^. It is easily seen that this is not an invariant concept, under modification of h 
in neghgible sets, while many applications of the theory to fluid dynamics (see for instance [23], 
|24) ) and conservation laws need this invariance property. This leads to the concept of regular 
Lagrangian flow (RLF in short): one may ask that, for all t £ [0,T], the image of the Lebesgue 
measure under the flow map x i— )• X{t,x) is still controlled by (see Definition 13. ip . 
It is not hard to show that, because of the additional regularity condition imposed on X, this 
concept is indeed invariant under modifications of b in Lebesgue negligible sets (see Remark 13. 8p . 
Hence RLF's are appropriate to deal with vector fields belonging to Lebesgue spaces. On the 
other hand, since this regularity condition involves all trajectories X{-,x) up to .^'^-negligible 
sets of initial data, the best we can hope for using this concept is existence and uniqueness of 
X{-,x) up to ^"'-negligible sets. Intuitively, this can be viewed as existence and uniqueness 
"with probability one" with respect to a reference measure on the space of initial data. Notice 
that already in the finite-dimensional theory different reference measures (e.g. Gaussian, see |5]) 
could be considered as well. 

To establish such existence and uniqueness, one uses that the concept of flow is directly 
linked, via the theory of characteristics, to the transport equation 

^/(s, x) + (6,(x), V,/(s, x)) = (18) 

and to the continuity equation 

^/xt + V • (btut) = 0. (19) 

The first equation has been exploited in jl^ to transfer well-posedness results from the transport 
equation to the ODE, getting uniqueness of RLF (with respect to Lebesgue measure) in W^. This 
is possible because the flow maps (s, x) i— )■ X(t, s, x) (here we made also explicit the dependence 
on the initial time s, previously set to 0) solve (jlSp for all t € [0, T]. In the present article, in 
analogy with the approach initiated in |T] (see also |16) for a stochastic counterpart of it, where 
(|19p becomes the forward Kolmogorov equation), we prefer rather to deal with the continuity 
equation, which seems to be more natural in a probabilistic framework. The link between the 
ODE (|17p and the continuity equation (|19p can be made precise as follows: any positive finite 
measure rj on initial values and paths, rj G x C([0,T];]R'^)), concentrated on solutions 

(x,7) to the ODE with initial condition x = 7(0), gives rise to a (distributional) solution to (|19p . 
with fit given by the marginals of r/ at time t: indeed, (|19p describes the evolution of a probability 
density under the action of the "velocity field" b. We shall call these measures rj generalized flows, 
see Definition 13.41 These facts lead to the existence, the uniqueness (up to ^'^-negligible sets) 
and the stability of the RLF X{t,x) in R'^ provided (O is well-posed in ([0, T]; L^nL°°(M'^)) . 
Roughly speaking, this should be thought of as a regularity assumption on b. See Remark 13.21 
and Section |6] for explicit conditions on b ensuring well-posedness. 

We shall extend all these results to flows on ^(M.'^^, the space of probability measures on 
M"^. The heuristic idea is that (|19p can be viewed as a (constant coefficients) ODE in the infinite- 
dimensional space ^{R'^), and that we can achieve uniqueness results for (|19p for "almost every" 
measure initial condition. We need, however, a suitable reference measure on ^(M'^), that we 
shall denote by u. Our theory works for many choices of v (in agreement with the fact that no 
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canonical choice of v seems to exist), provided v satisfies the regularity condition 

see Definition 13.51 (See also Example 13.61 for some natural examples of regular measures u.) 
Given u as reference measure, and assuming that (|19p is well-posed in L'^ ([0,T]; nL°°(M'^)), 
we prove existence, uniqueness (up to i/'-negligible sets) and stability of the regular Lagrangian 
flow of measures fx. Since this assumption is precisely the one needed to have existence and 
uniqueness of the RLF X{t,x) in R'^, it turns out that the RLF fj,{t,fi) in J^(W^^ is given by 

Kt,fi)= [ Sxit,.) dfii^) V t G [0, T] , ^ G ^ {W") , (20) 

which makes the existence part of our results rather easy whenever an underlying flow X in 
M"^ exists. On the other hand, even in this situation, it turns out that uniqueness and stability 
results are much stronger when stated at the =f^(M'^) level. 

In our proofs, which follow by an infinite-dimensional adaptation of [T], we use also the 
concept of generalized flow in ^(W^^ , i.e. measures r] on ^{R'^) xC([0,T]; ^(R'^)) concentrated 
on initial data/solution pairs {fi,uj) to (|19p with uj{0) = fJ., see Definition 13.91 

Acknowledgement. We thank Dr. Marilena Ligabo for pointing out a serious gap in a pre- 
liminary version of Theorem 17.11 The first author was partially supported by ERG ADG Grant 
GeMeThNES and the second author was partially supported by NSF Grant DMS-0969962. 



2 Notation and preliminary results 

Let X be a Polish space (i.e. a separable topological space whose topology is induced by a 
complete distance). We shall denote by B{X) the cj-algebra of Borel sets of X, by ^(X) (resp. 
^(X), ^+(X)) the space of Borel probability (resp. finite Borel, finite Borel nonnegative) 
measures on X. For A G B{X) and u G ^(X), we denote by ul^A G ^(X) the restricted 
measure, namely iyLA{B) = u{A n B). Given f : X ^ Y Borel and /i G ^(X), we denote 
by /jA* ^ ^[Y^ the push-forward measure on Y, i.e. /jj//(A) = fi{f~^{A)) (if ^ is a probability 
measure, f^/i is the law of / under fi) and we recall the basic integration rule 

/ (j) df^fj, = (/) o f dfi 4> bounded and Borel. 
Jy Jx 

We denote by XA the characteristic function of a set A, equal to 1 on A, and equal to on its 
complement. Balls in Euclidean spaces will be denoted by Bji^xq), and by Bji if xq = 0. 

We shall endow ^(X) with the metrizable topology induced by the duality with Cfe(X), the 
space of continuous bounded functions on X: this makes =^^(X) itself a Polish space (see for 
instance |31 Remark 5.1.1]), and we shall also consider measures v G ^+(=f^(X)). 

Typically we shall use greek letters to denote measures, boldface greek letters to denote 
measures on the space of measures, and we occasionally use d^ for a bounded distance in J^(X) 
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inducing the weak topology induced by the duahty with Cb{X) (no specific choice of will be 
relevant for us). We recall that weak convergence of /i„ to ^ implies 



lim / f dfin = f dfi for all / bounded Borel, with a /i- negligible discontinuity set. 
Jx Jx 

(21) 

Also, in the case X = M'^, recall that a sequence (/x„) C ^{M.'^^ weakly converges to a probability 
measure fi in the duality with Cfe(M'^) if and only if it converges in the duality with (a dense 
subspace of) Cc(M'^). 

We shall consider the space C([0,T]; ,^(M'^)), whose generic element will be denoted by u!, 
endowed with the sup norm; for this space we use the compact notation ^^(^(M'^)). We also 
use Ct as a notation for the evaluation map at time t, so that et{u}) = uj{t). Again, we shall 
consider measures r] G ^+(il7'(^(R'^))) and the basic criterion we shall use is the following: 

Proposition 2.1 (Tightness). Let (r/„) C ^-^.[Qt{^0^'^))) be a bounded family satisfying: 



(i) (space tightness) for all e > 0, sup„ riA{^uj: sup ti;(i)(IR'^ \ Br) > e} j — )• as — >• oo; 

(ii) (time tightness) for all (j) G C'^{W^), n > 1, the map 1 1— )• j^^i (/^duj^t) is absolutely continu- 
ous in [0,T] for Tjj^-a.e. uj and 



lim sup r]„[ \uj 

Mtoo n 



(j) duj{t) 



' dt> M 



})-■ 



Then {rj^) is tight. 

Proof. For all (/> G C^{W^) we shall denote by : 9.t{^{W^)) C([0,r]) the time-dependent 
integral w.r.t. (j). Since the sets 

|/GH^i'i(0,r): sup|/|<C, £\f'{t)\dt<M 

are compact in C([0,r]), by assumption (ii) the sequence {{I(f,)^T]n) is tight in ^+(C([0,T])) 
for all (j) S C^{W^). Hence, if we fix a countable dense set (0^) C C^(M'^) and e > 0, we can 
find for /c > 1 compact sets iTf C C{[0,T]) such that sup„ t/„(Qt(^(M'^)) \ I^j^i^k)) < e2-^ . 
Thus, if denotes the intersection of all sets I^^{Kf,), we get 



supr/„(Or(^(M'^))\K^)<e. 

n 

Analogously, using assumption (i) we can build another compact set L'' C Ot(=5^(M'^)) such 
that sup„ T7„(i7T(<^(]R'^)) \ i^) < e and, for all integers A: > 1, there exists R = such that 
oj{t){W^ \ Br) < 1/k for all uj e and t e [0,T]. 

In order to conclude, it suffices to show that K"^ n is compact in if ('^p) C 

K"^ n we can use the inclusion in I^^{Kf,) and a diagonal argument to extract a subsequence 
{^p{e)) such that J (j)^ dujp(^£'j{t) has a limit for all t G [0,T] and all /c > 1 and the limit is 
continuous in time. By the space tightness given by the inclusion (ujp) C L^, ujp(^£'j{t) converges 
to uj{t) in ^(M'^) for all t G [0, T], and t ^ u;{t) is continuous. □ 
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The next lemma is a refinement of [21 Lemma 22] and |30| Corollary 5.23], and allows to 
obtain convergence in probability from weak convergence of the measures induced on the graphs. 

Lemma 2.2. Let fn '■ X ^ Y , f : X ^ Y be Borel maps, Un, v G !^{X) and assume 
that (Id X fn)ii'n weakly converge to (Id x in X x Y . Assume in addition that we have 

the Skorokhod representations m„ = (i„)jP, = i^P, with {W,J-',¥) probability measure space, 
in, i '■ W ^ X measurable, and in — )• i ¥-almost everywhere. 
Then fn°in^f°i if^ "^-probability. 

Proof. Let dy denote the distance in Y. Up to replacing dy by minjdy,!}, with no loss of 
generality we can assume that the distance in Y does not exceed 1. Fix e > and g £ Cb{X; Y) 
with J-^ dy{g, f) du < e^. We have that {dy{fn ° in, f ° i) > 3e} is contained in 

{dy{fn °in,g° in) > e} U {dy{g o in,g o i) > e}U {dy{g oi,foi)> e}. 

The second set has infinitesimal P-probability, since g is continuous and i„ — )■ i P-a.e.; the third 
set, by Markov inequality, has P-probability less than e; to estimate the P-probability of the first 
set we notice that 

f{{dy{fn o in, go in) > e}) = ^n{{dy{fn,g) >£})<-[ xd{Id X /„)jjl/„ 

with xi^,y) '■= dy{g{x),y). The weak convergence of (Id x fn)^Vn yields 

limsupP({(iy(/„ o o i„) > e}) < - / xd{Idxf)>^v 

= - I dy{g{x),f{x))du{x)<e. 
£ Jx 

□ 

3 Continuity equations and flows 

In this section we shall specify the basic assumptions on b used throughout this paper, and the 
conventions about (|19p concerning locally bounded respectively measure-valued solutions. We 
shall also collect the basic definitions of regular flows we shall work with, recalling first those 
used when the state space is and then extending these concepts to ^(M'^). 

3.1 Continuity equations 

We consider a Borel vector field b : [0,T] x M.'^ —?■ M.'^, and set fet(-) := b{t, ■); we shall not work 
with the Lebesgue equivalence class of b, although a posteriori our theory is independent of the 
choice of the representative (see Remark 13. 8p : this is important in view of the fact that (|19p 
involves possibly singular measures. Also, we shall not make any integrability assumption on 
b besides LI^^(^[0,T] x M.'^) (namely, the Lebesgue integral of |6| is finite on [0,T] x Bji for all 
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R > 0); the latter is needed in order to give a distributional sense to the functional version of 
(|19p . namely 

j^wt + V- (btwt) = (22) 

coupled with an initial condition wq = w ^ L^^{M.'^), when wt is locally bounded in space-time. 

It is well-known and easy to check that any distributional solution w^tjx) = Wt{x) to (|22p 
with Wt locally bounded in uniformly in time, can be modified in a ^^-negligible set of times 
in such a way that t ^ wt is continuous w.r.t. the duality with Cc(lR'^), and well-defined limits 
exist at t = 0, t = T (see for instance [H Lemma 8.1.2] for a detailed proof). In particular 
the initial condition wq = w \s then well defined, and we shall always work with this weakly 
continuous representative. 

In the sequel, we shall say that the continuity equation (|22p has uniqueness in the cone of 
functions ([0, T]; n L°°(M'^)) if, for any w £ L^ClL'^iR'^) nonnegative, there exists at most 
one nonnegative solution wt to (|22p in L°°([0,T];L^ n L°°(M'^)) satisfying the condition 

uiQ = w. (23) 

Coming to measure- valued solutions to (|19p . we say that t € [0, T] i— ?■ //t G ^+(M'^) solves (|19p 
if |6| € Ly^^(^(0,T) X Mf^ ; fifdt) , the equation holds in the sense of distributions and t^f cl)d^t 
is continuous in [0, T] for all ^ G Cc(M'^). 

3.2 Flows in 

Definition 3.1 (z^-RLF in W^). Let X : [0,r] xW^ and v G ^+(M'^) with v ^'^ and 

with bounded density. We say that X is a v-RLF in (relative to b £ -^ioc((^' ^ ■'^'^) ) ^/ 
following two conditions are fulfilled: 

(i) for p-a.e. x, the function t i— )■ X{t,x) is an absolutely continuous integral solution to the 
ODE dm in [0,T] with X{0,x) = x; 

(ii) X{t, ■)^u < C^^ for all t G [0,T], for some constant C independent oft. 

Notice that, in view of condition (ii), the assumption of bounded density of v is necessary 
for the existence of the z^-RLF, as X(0, •)jji/ = v. 

In this context, since all admissible initial measures v are bounded above by C££'^ ^ uniqueness 
of the 1/-RLF can and will be understood in the following stronger sense: ii f, g £ L^{W^) Pi 
L°°(M'^) are nonnegative and X and Y are respectively a /^"^-RLF and a gJ("^-BIjF, then 
X{-, x) = Y{-,x) for ^'^-a.e. x £ {f > 0} n {g > 0}. 

Remark 3.2 {BV vector fields). We shall use in particular the fact that the I'-RLF exists for all 
u < C^"^, and is unique, in the strong sense described above, under the following assumptions on 
b: \b\ is uniformly bounded, bt G 5yioc(M'^; M'') and V ■ bt = gt^'^ < for if^-a.e. t G (0,r), 
with 

||5t||Loo(Rd) G Li(0,r), \Dbt\{BR) £ L\Q,T) for all ii > 0, 

where \Dbt\ denotes the total variation of the distributional derivative of bt. See |T] or |2] and 
the paper |12] for Hamiltonian vector fields. 
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Remark 3.3 (^'^-RLF). In all situations where the z^-RLF exists and is unique, one can also 
define by an exhaustion procedure a J^'^-KLF X, uniquely determined (and well defined) by the 
property 

X{-, x) = Xf (•, x) ^'^-a.e. on {/ > 0} 

for all / G L°° n L^(R'^) nonnegative, where X^ is the /^'^-flow. Also, it turns out that if (|22]) 
has backward uniqueness, and if the constant C in Definition 13. If ii) can be chosen independently 
of < then X{t, ■)^^'^ < C5£'^ . We don't prove this last statement here, since it will not 
be needed in the rest of the paper, and we mention this just for completeness. 

In the proof of stability and uniqueness results it is actually more convenient to consider a 
generalized concept of flow, see |2] for a more complete discussion. We denote the evaluation 
map (x,uS) G M'^ X C([0,r];M'^) ^ uj{t) G R'^ again with ej. 

Definition 3.4 (Generalized z^-RLF in M''). Let v G ^+{W^) and rj G ^(M'^ x C([0, T]; M'^)) . 
We say that r] is a generalized u-RLF in (relative to h) if: 

(i) ieo)tV = ^; 

(a) T] is concentrated on the set of pairs {x,j), with 7 absolutely continuous solution to (|17p . 
and 7(0) = x; 

(Hi) {et)^r] < C^'^ for all t G [0, T], for some constant C independent oft. 
3.3 Flows in ^{W^) 

Given a nonnegative cr- finite measure u G ^+ (=f^(R'^)) , we denote by Ei/ G ^+(M'^) its expec- 
tation, namely 

/ (j)d¥.u = I I (()d^du{ii) for all bounded Borel. 

Definition 3.5 (Regular measures on ^+(^(M'^))). Let u G ^+(^(M'^)). We say that v is 
regular ifW,v < CJi"^ for some constant C. 

Example 3.6. (1) The first standard example of a regular measure u is the law under p^'^ of 
the map x 6x, with p G L^(M.'^) n L°°{W^) nonnegative. Actually, one can even consider the 
law under and in this case u would be cr-finite instead of a finite nonnegative measure. 

(2) If d = 2n and z = {x,p) G x M" (this factorization corresponds for instance to 
flows in a phase space), one may consider the law under p^'^ of the map a; 1— )• x 7, with 
p G L^{W^) n L°°(]R") nonnegative and 7 G ^{Rp) with 7 < C^"; one can also choose 7 
dependent on x, provided x 1— >• 7^; is measurable and < for some constant C independent 
of X. 

(3) We also expect that the entropic measures built in [28], |29) are regular, see also the 
references therein for more examples of "natural" reference measures on the space of measures. 
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As we explained in the introduction, Definition 13.11 has a natural (but not perfect) transpo- 
sition to fiows in {Mf) : 

Definition 3.7 (i/-RLF in ^(M'^)). Let fx : [0,T] x ^{R'^) ^{R'^) and v G ^+(=^(R'^)). 
We say that n is a u-RLF in ^{W^) (relative to b with \b\ £ Lj^Xi^^ ^) ^ ^'^^ t-kdt)) if 

(i) foru-a.e. fi,t^ fXt '■= /^(^i/^) (weakly) continuous from [0,T] to ^3^(R'^^ with /i(0,/i) = 
/i and fit solves (|19p in the sense of distributions; 

(ii) E{fi{t, Op) < C^'^ for all t G [0,r], for some constant C independent oft. 

Notice that no I'-RLF can exist if i/ is not regular, as yLt(0, Op = Notice also that 
condition (ii) is in some sense weaker than /^(t, Op < Cv (which would be the analogue of (ii) in 
Definition 13. II if we were allowed to choose u = see also Remark 13. 3p . but it is sufficient for 
our purposes. As a matter of fact, because of infinite-dimensionality, the requirement of quasi- 
invariance of u under the action of the flow /x (namely the condition fj,{t, Op ^ i^) would be a 
quite strong condition: for instance, if the state space is a separable Banach space V, the reference 
measure 7 is a nondegenerate Gaussian measure, and b{t,x) = v, then X{t,x) = x + tv, and 
the quasi-invariance occurs only if v belongs to the Cameron-Martin subspace H of V, a dense 
but 7-negligible subspace. In our framework. Example I3.6f 2) provides a natural measure u that 
is not invariant, because its support is not invariant, under the flow: to realize that invariance 
may fail, it suffices to choose autonomous vector fields of the form b{x,p) := (p, — VC/(x)). 

Remark 3.8 (Invariance of v-KLF). Assume that fj,{t,fi) is a 1/-RLF relative to b and 6 is a 
modification of b, i.e., for ^^-a.e. t G (0,T) the set Nt := {bt 7^ bt} is ^'^-negligible. Then, 
because of condition (ii) we know that, for all t G (0, T), fx{t, fi){Nt) = for iz-a.e. fi. By 
Fubini's theorem, we obtain that, for u-a.e. fi, the set of times t such that fx{t, fi){Nt) > is 
^^-negligible in (0, T). As a consequence 1 1— )■ /x(t, /u) is a solution to (|19p with bt in place of bt, 
and fj. is a t^'-RLF relative to b as well. 

In the next definition, as in Definition 13. 4[ we are going to consider measures on ^(R'^) x 
Ot(,^3^(M'^)), the first factor being a convenient label for the initial position of the path (an 
equivalent description could be given using just measures on Qri^O^'^)), at the price of an 
heavier use of conditional probabilities, see |2l Remark 11] for a more precise discussion). We 
keep using the notation et for the evaluation map, so that et{fJ.,uj) = Lj{t). 

Definition 3.9 (Generalized u-RLF in ^{R'^)). Let v G ^+(^(M'^)) and r] G ^+{^{R'^) x 
Qt{^{^'^))) ■ We say that rj is a generalized u-RLF in ^^(R'^^ (relative to b with \b\ G 
Ll^{{0,T)xR'';fitdt)) if: 

(i) (eo)tt^7 = v; 

(ii) T] is concentrated on the set of pairs {fi,uj), with u solving (|19p . uj{0) = fi.; 
(ill) E((et)tt77) < C^'^ for all t G [0,r], for some constant C independent oft. 
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Again, by conditions (i) and (iii), no generalized zv-RLF can exist if v is not regular. Of 
course any I/-RLF yii induces a generalized i/-RLF r}: it suffices to define 

V ■■= (24) 

where 

: ^(M'^) ^ ^(R'^) X Qri^iR'')), *^(^) := (^i, /x)). (25) 

It turns out that existence results are stronger at the RLE level, while results concerning unique- 
ness are stronger at the generalized RLF level. 

The transfer mechanisms between generalized and classical flows, and between flows in =f^(R'^) 
and flows in M.'^ are illustrated by the next proposition. 

Proposition 3.10. Let r] be a generalized u-RLF in ^(R'^) relative to b. Then: 
(i) Er/ is a generalized Ku-RLF in R'^ relative to b; 
(a) the measures fit '■= IE((et)jr/) = (ej)jjE77 G ^+(R'^) satisfy (fT9]) . 
In addition, fit = Wt^'^ with w E L'^ {[0,T]; L^ D L°°(R'^)). 

Proof. Statement (i) is easy to prove, since the continuity equation is linear. Statement (ii), 
namely that (single) time marginals of generalized flows in R*^ solve (|19p , is proved in detail in 
[21 Page 8]. The final statement follows by the regularity condition on r/. □ 

4 Existence and uniqueness of regular Lagrangian flows 

In this section we recall the main existence and uniqueness results of the i^-RLF in R"^, and see 
their extensions to i^-RLF in ^(R'^). It turns out that existence and uniqueness of solutions 
to (j22p in L5^([0, T];L^ nL°°(R'^)) yields existence and uniqueness of the v-RLF, and existence 
of this flow implies existence of the iz-RLF when iv is regular. Also, the (apparently stronger) 
uniqueness of the i^-RLF is still implied by the uniqueness of solutions to (|22p in L'^ ([0, T];L^n 
L°°(R'^)). 

The following result is proved in [21 Theorem 19] for the part concerning existence and in [21 
Theorem 16, Remark 17] for the part concerning uniqueness. 

Theorem 4.1 (Existence and uniqueness of the z^-RLF in R*^). Assume that (|22p has existence 
and uniqueness in L^ ( [0, T];L^n L°° (R*^)) . Then, for all v G ^{W^) with u <^ ^'^ and bounded 
density the v-RLF in exists and is unique. 

Now we can easily show that existence of the z^-RLF implies existence of the 1/-RLF, by a 
superposition principle. However, one might speculate that, for very rough vector fields, a 1/-RLF 
might exist in ^(R'^), not induced by any z^-RLF in R"^. 
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Theorem 4.2 (Existence of the v-RLF in ^(M*^)). Let v e Ji{W^) with v < 5^'^ and bounded 
density, and assume that a v-RLF X in exists. Then, for all v S ^^{^S^iW^y) with Ku = v, 
a v-RLF /X in ^(M'^) exists, and it is given by 

fj,{t,fi):= Sx(t,x)dKx). (26) 

Proof. The first part of property (i) in Definition 13.71 is obviously satisfied, since the fact that 
t I— 7- X{t,x) solves the ODE for some x corresponds to the fact that t i— t- Sx{t,x) solves (|19p . 
On the other hand, since iv is regular and X is a RLE, we know that X{-,x) solves the ODE 
for Ei>'-a.e. x; it follows that, for u-a.e. n, X(-,x) solves the ODE for ^u-almost every x, hence 
H{t,fi) solves (|19p for v-a.e. fi. This proves (i). 
Property (ii) follows by 

/ (f){x)dE{n{t,-)f^u){x) = / 4)diJi{t,^)du{p) 

(j){X{t, x)) dfi{x) dv>{fi) 



= [ cp{X{t,x))dv{x) <CL j (P{z)dz 

JRd jRd- 

where C is the same constant in Definition I3.1fii) and L satisfies u < LJ^f^. □ 



The following lemma (a slight refinement of [U Theorem 5.1] and of [5l Lemma 4.6]) provides a 
simple characterization of Dirac masses for measures on Cw{[0, T];E^ and for families of measures 
on E. Here ii^ is a closed, convex and bounded subset of the dual of a separable Banach space, 
endowed with a distance ds inducing the weak* topology, so that (E,dE) is a compact metric 
space; C^([0,T];£^) denotes the space of continuous maps with values in {E,dE), endowed with 
sup norm (so that these maps are continuous with respect to the weak* topology). We shall 
apply this result in the proof of Theorem 14.41 with 

E := {/i G ^{R'^) : I^KM"^) < l} I) ^{R'^), (27) 

thought as a subset of (Co(M'^))*, where Co(M'^) denotes the set of continuous functions vanishing 
at infinity (i.e. the closure of Cc(M'^) with respect to the uniform convergence). 

Lemma 4.3. Let E C G* , with G separable Banach space, be closed, convex and bounded, and 
let a be a positive finite measure on Ct„ ([0, T]; ii^) . Then a is a Dirac mass if and only if {et)^(T 
is a Dirac mass for all t £ QCi [0, T] . 

If {F,T ,\) is a measure space, and a Borel family {i'z}z(^F of probability measures on E (i.e. 
-2 ^ T^z{A) is T -measurable in F for all A d E Borel) is given, then are Dirac masses for 
X-a.e. z £ F if and only if for all y £ G and c G M there holds 

u^{{x£E: {x,y) < c})u^{{x £ E : {x,y)>c}) = for X-a.e. z £ F. (28) 
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Proof. The first statement is a direct consequence of tlie fact that all elements of Cu,{[0,T]; 
are weakly* continuous maps, which are uniquely determined on Q H [0,T]. In order to prove 
the second statement, let us consider the sets Ay := {x £ E : {x,yi) < Cj}, where yi vary in a 
countable dense set of G and cj are an enumeration of the rational numbers. By (|28p we obtain 
a A-negligible set Nij £ satisfying iyz{Aij)i'z{E\Aij) = for all z G F\Nij. As a consequence, 
each measure z^^, as z varies in F \ Nij, is either concentrated on Aij or on its complement. 
For z £ F \ UjNij it follows that the function x i— )• {x,yi) is equivalent to a constant, up to 
i/^-negligible sets. Since the functions x i— )■ {x, yi) separate points of E, is a Dirac mass for all 
z £ F \ UijNij as desired. □ 

The next result shows that uniqueness of (HH]) in L^{[0,T];L^ n L°°(R'^)) and existence of 
a generalized 1/-RLF imply existence of the 1/-RLF and uniqueness of both, the I/-RLF and the 
generalized i^-RLF. 

Theorem 4.4 (Existence and uniqueness of the u-RLF in Assume that (p2|) has 

uniqueness m L^([0,r];Li n L°°(M'^)). // a generalized u-RLF in g^{W^) rj exists, then the 
v-RLF fi in ^(M'^) exists. Moreover they are both unique, and related as in (|24p . (|25p . 

Proof. We fix a generalized i^-RLF rj and we show first that r/ is induced by a I/-RLF (this 
will prove in particular the existence of the iv-RLF). To this end, denoting by vr : ^(R'^) x 
Qt{^{^'^)) — ^ ^(W^) the projection on the first factor, we define by 

rj^ ■.= E{ri\7T = fi)£^{nT{^{R% 

the induced conditional probabilities, so that dr]{fi,uj) = dr]^{uj)diy{fi). Taking into account the 
first statement in Lemma |4.3[ it suffices to show that, for t £ Qf] [0,T] fixed, the measures 

9^ := E((e,-)jT7|o^(0) = /i) = (e,-)jr,^ £ ^+(^(M'^)) 

are Dirac masses for u-a.e. jj, £ ^(R'^). Still using Lemma [4. 3 1 we will check the validity of (|28p 
with \ = v. Since 0^ = 5^ when t = 0, we shall assume that t > 0. 

Let us argue by contradiction, assuming the existence of L G B{^{R'^)) with u{L) > 0, 
£ Co(R'^), c G R such that both e^{A) and 6'^(^(R'^) \ A) are strictly positive for all n £ L, 
with 

A:= !^p£ ^{M.'^) : j 0dp<c|. 

We will get a contradiction with the assumption that the equation (|22p is well-posed in ([0, T]; L^n 
L°°(R'^)), building two distinct nonnegative solutions of the continuity equation with the same 
initial condition w £ r\ L°°(R'^). With no loss of generality, possibly passing to a smaller set 
L still with positive i/-measure, we can assume that the quotient ^(/i) := ^{A) / 9 ^(^3^ (W^^ \ ^4) 
is uniformly bounded in L. Let C ^^^(^(R'^)) be the set of trajectories oj which belong to A 
at time i, and let Q.2 be its complement; we can define positive finite measures ry*, i = 1, 2, in 
^(R'') X ^t{^{^'^)) by 

d-n^{ix,uj) := d{x^^'n^){u:)d{xLJ^){lj), drf{ii,u:) := d{xn2Vf^){(^)d{xL9i^){lJ')- 
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By Proposition 13.101 both rj^ and r]^ induce solutions vj-j:, to the continuity equation which are 
uniformly bounded (just by comparison with the one induced by 77) in space and time. Moreover, 
since 

and analogously 

ieo)ir,^ = e,{^{R^)\A)xL{f^)g{fi)iy, 

our definition of g gives that (60)^^7"^ = {eo)^'^. Hence, both solutions w^, w1 start from the 
same initial condition iv{x)^ namely the density of E(0^(A)xl(^)i^) with respect to On the 
other hand, it turns out that 



(t)w\dx =111 (j)duj{t) dr]^{io) duip) 

XAi^it)) / (t)duj{t)drf^{uj)du{fi) 




(t)dpde^{p)du{^i) <c / 9f,{A)di^{fi). 

Analogously, we have 

(t>w} dx>c I e^{^{R'^) \ A)g{f,) du{ii) = c I e^{A) du{p). 
Jl Jl 

Therefore 7^ and uniqueness of the continuity equation is violated. 

Now we can prove uniqueness: if a is any other generalized i^-RLF, we know a is induced 
by a 1/-RLF, hence for iv-a.e. p also the measures E(<t|u;(0) = //) are Dirac masses; but, since 
the property of being a generalized flow is stable under convex combinations, also the measures 
(corresponding to the generalized iv-RLF (ry + <t)/2) 

^E(77|a;(0) = p) + ^E(cr|a;(0) = = E (^^?-±-^|a;(0) = ^ 

must be Dirac masses for i/-a.e. ^. This can happen only if E(77|a;(0) = p) = E(<t|w(0) = p) 
for f-a.e. /i, hence a = rj. Finally, since distinct v-KLF fi and fi' induce distinct generalized 
f-RLF r] and r]', uniqueness is proved also for i^-RLF. □ 



5 Stability of the u-RLF in ^{R'^) 

In the statement of the stability result we shall consider varying measures G ^(^(M'^)), 
n > 1, and a limit measure v. (The assumption that all i>n are probability measures is made in 
order to avoid technicalities which would obscure the main ideas behind our stability result, and 
one can always reduce to this case by renormalizing the measures. Moreover, in the applications 
we have in mind, our measures i>n will always have unitary total mass.) We shall assume that the 
Vn are generated as {in)(,R, where (VF, P) is a probability measure space and in '■ W ^ ^(W^^ 
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are measurable; accordingly, we shall also assume that v = ijjP, with in ^ i P-almost everywhere. 
These assumptions are satisfied in the applications we have in mind, and in any case Skorokhod's 
theorem (see §8.5, Vol. II]) could be used to show that weak convergence of to v always 
implies this sort of representation, even with W = [0, 1] endowed with the standard measure 
structure, for suitable in, i- 

Many formulations of the stability result are indeed possible and we have chosen one specific 
for the application we have in mind. Henceforth we fix an autonomous vector field & : M"' — )• 
satisfying the following regularity conditions: 

(a) d = 2n and b{x,p) = {p, c{x)), {x,p) G M'^, c : M" — )• M" Borel and locally integrable; 

(b) there exists a closed ^"-negligible set S such that c is locally bounded on M" \ 5; 

(c) the discontinuity set S of c is ^"-negligible. 

Lemma 5.1. Let S C M" closed, and assume that b is representable as in (a) above. Let 
fit : [0,T] —7- ^(M'^) be solving (|19p in the sense of distributions in (M" \ S) x M" and assume 
that ^ 

/ / 1-78} — 7:7 dnt{x,p)dt <oo Vii>0 
Jo J Br disf^^x, b) 

for some (3 > 1 (with the convention 1/0 = +oo^. Then (|19p holds in the sense of distributions 
in W^. 

Proof. First of all, the assumption implies that x M") = for ^^-a.e. t G (0, T). The proof 
of the global validity of the continuity equation uses the classical argument of removing the 
singularity by multiplying any test function (j) £ C^i^W^) by Xk, where Xk{x) = x(fcdist(x, /S)) 
and X is a smooth cut-off function equal to on [0, 1] and equal to 1 on [2, -|-oo), with < x' ^ 2. 
If we use (/)Xfc as a test function, since Xk depends on x only, we can use the particular structure 
(a) of b to write the term depending on the derivatives of Xk as 

^ / / <^x'(^dist(x, S'))(p, Vdist(x, 5)) (i;Uj (x,p)(it. 

If K is the support of (j), the integral above can be bounded by 

, n f , , / N , 2^^^ max/^ \p(/)\ f 1 , , ^ , 

2max|p(/)| / / kdfit[x,p)dt < j / / „ — — dfj.t[x,p)dt 



^ Jo J{x£K:kdist{x,S)<2} ' ~ Jo ii^ dist'^ (x, 5) 

and as ^ > 1 the right hand side is infinitesimal as A; — )■ oo. □ 

The following stability result is adapted to the application we have in mind: we shall apply 
it to the case when //^(t, /i) are Husimi transforms of wavefunctions. 

Theorem 5.2 (Stability of the v-RLF in ^(R'^)). Let in, i be as above and let fin '■ [0;^] ^ 
in{W) — )■ ^(W^) be satisfying /^„(0, i„(t(;)) = in{w) and the following conditions: 
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(i) (asymptotic regularity) 



lim sup 

n—^oo J W 



cl)dn^{t,in{w))dF{w) <C (pdx 



for all (f) G Cc(M ) nonnegative, for some constant C independent of t; 
(ii) (uniform decay away from the singularity) for some /3 > 1 

1 



sup lim sup 

5>0 n— >oo 



W Jo JBa dist^(x, 5) + (5 



dfj.j^{t,in{w)) dt dF{w) < oo 



\/R > 0; 



(29) 



as 



(Hi) (space tightness) for all 6 > 0, F({w e W : sup fi^{t,in{w)){M.^ \ Br) > d]) 

^ te[o,T] ' 

i? — 7- oo uniformly in n; 

(iv) (time tightness) for ¥-a.e. w ^ W , for all n > 1 and (j) ^ {W^) , t ^ f^a dfJ-nitjiniw)) 
is absolutely continuous in [0, T] and, uniformly in n, 



'{{' 



lim P( <! u; G 1^ 

Affoo 

(v) (limit continuity equation) 



T 



(j)dtJ,n{t,in{w)) 



dt > M 



}) 



lim 

n— >oo 



W 



ip'{t) I (f>dnn{t,in{w)) + ip{t) (5, V(/>) 



dt 



dF{w) = 

(30) 



for all (f> G C7^^(IR'^\ (5 x M")), if G C^{0,T). 

Assume, besides (a), (b), (c) above, that (|22p has uniqueness m ([0, T]; n L°°(M'^)) . Then 
the v-RLF fj,{t, fi) relative to b exists, is unique, and 



lim 

n— >oo 



sup d^{Hn{t,in{w)), fi{t,i{w))) dF{w) = 0. 

iyt6[0,T] 



(31) 



Proof Let {rjj C ^+(^(M'^) x ^Ti^O^ ))) be induced by /i„ pushing forward iv„ = (i„)ftP 
via the map i— )• (/x, ^„(t, /i)). Conditions (iii) and (iv) correspond, respectively, to conditions 
(i) and (ii) of Proposition 12. 1| hence the marginals of r)^ on Qt(c5^(R'^)) are tight; since the first 
marginals, namely Un, are tight as well, a simple tightness criterion in product spaces (see for 
instance |H Lemma 5.2.2]) gives that (77^) is tight. We consider a weak limit point ij of (77^) and 
prove that r] is the unique generalized i^-RLF relative to b, this will give that the whole sequence 
iVn) weakly converges to t]. Just to simplify notation, we assume that the whole sequence (77^) 
weakly converges to 77. 

We check conditions (i), (ii), (iii) of Definition 13.91 First, since /x„(0,^) = /.f i/„-a.e., we get 
(eo)j|'/7„ = I'n, hence {eo)^ = v and condition (i) is satisfied. Second, we check condition (iii): 
for (j) G Cc(M'^) nonnegative we have 



())dE{{et)ir]n 



)(iE(/i(t, •)jji/^ 



w 



<j)dfj,n{t,iniw))dF{w) 
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and we can use assumption (i) to conclude that 



(32) 



so that condition (iii) is fulfilled. 

Finally we check condition (ii). Since r/„ are concentrated on the closed set of pairs (/U,a;) 
with 1^(0) = /.i, the same is true for 77; it remains to show that uj{t) solves (|19p for T^-a.e. (fj,,uj). 
We shall denote by <t G ^^(^Qt{^{^'^))) the projection of rj on the second factor and prove 
that pop holds for cr-a.e. uj. 

We fix (/) G (M'^ \ (5 X R")) and (f G 0^(0, T); we claim that the discontinuity set of the 
bounded map 



is <T-negligible. Indeed, using (|2ip with X 
contained in 



(p'{t) / (f)du}{t) + ip{t) / {b,V(j)) du}{t) 



dt 



(33) 



this discontinuity set is easily seen to be 



u E Ot(^(M'')) : ^ a;(t)(S x W) > o| , 



(34) 



where S is the discontinuity set of c. Since ^'^{Ti x M") = 0, by assumption (c), for all t € [0, T] 
the inequality ([32]) gives tj(t)(SxM") = for cr-a.e. w; by Fubini's theorem in [{),T]xQ.t{^{'^'^)) 
we obtain that the set in (j34)) is cr-negligible. 

Now we write assumption (|29p in terms of r/„ as 



sup lim sup 

5>0 n-5.oo J3^(Rd)xCT(i32{K'*)) -^0 



1 



■ doj{t) dtdr]n{fj.,uj) < 00 Vi? > 0, 



iBa dist'^(a;,5) + 5 

and take the limit thanks to Fatou's Lemma and the Monotone Convergence Theorem to obtain 

1 



nT(.i^{Rd)) Jo J Br dist''^(x, S) 
Next we write assumption (v) in terms of rj^ as 



■ du {t)dt do- {u}) < 00 Vi?>0. 



(35) 



lim / 



C 



ip'{t) / (j)doo{t) + ip{t) / {b,V(l))duj{t) 



dt 



C?^n(/^>W) = 



with C E Cf,(^(]R"') X Qt{^{^'^))) nonnegative; then, the claim on the continuity of the map 
in dMI) and ([21$ with X = ^{W^) xQt{^(^'^)) allow to conclude that 



c 



ip'{t) I 4>dio{t) + ip{t) I {h,V4))duj{t) 



dt 



dr]{fi,uj) = 0. 



Now we fix ^ C C^lW^ \ (5 x M")), B C C^{0,T) countable dense, and use the fact that C is 
arbitrary to find a <T-neglig ible set N C 17t(^(M'^)) such that 



T 



if' it) I 4>doj{t) + ip{t) I {bt,V(p) dtj{t) 



dt = o \f(j)eA,'iipeB 
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for all uj ^ N, and by a density argument we conclude that a is concentrated on solutions to 
the continuity equation in M'^ \ (5 x W^). By Lemma l5.ll and (|35p we obtain that cr-a.e. the 
continuity equation holds globally. 

By Theorem 14.41 we know that the u-RLF in ^(W^^ exists, is unique, and related to 

the unique generalized iv-RLF rj as in (|24p . (|25p . This proves that we have convergence of the 
whole sequence (t7„) to rj. By applying Lemma |2^ with X = J^(M.'^) and Y = ^Ti^O^"^)) we 
conclude that (|3ip holds. □ 

In the next remark we consider some extensions of this result to the case when 6 satisfies (a) , 
(b) only, so that no information is available on the discontinuity set S of c. 

Remark 5.3. Assume that b satisfies (a), (b) only. Then the conclusion of Theorem 15.21 is still 
valid, provided the asymptotic regularity condition (i) holds in a stronger form, namely 

/ / (t)dfi^{t,in{w))dF{w) < C ct)dx V G Cc(M'^), > 0, n > 1 

for some constant C independent of t. Indeed, assumption (c) was needed only to pass to the 
limit, in the weak convergence of to r/, with test functions of the form (|33p . But, if the 
stronger regularity condition above holds, convergence always holds by a density argument: first 
one checks this with b continuous and bounded on suppi;^, and in this case the test function 
is continuous and bounded; then one approximates b in on suppc/) by bounded continuous 
functions. 



6 Well-posedness of the continuity equation with a singular po- 
tential 

In this section we shall assume that d = 2n and consider a more particular class of autonomous 
and Hamiltonian vector fields 5 : M*^ — )• R"' of the form 

b{z) = {p,-VU{x)), z = G X M". 

Having in mind the application to the convergence of the Wigner/Husimi transforms in quantum 
molecular dynamics, we assume that: 

(i) there exists a closed ^"-negligible set S C M" such that U is locally Lipschitz in M" \ S 
and VU G SMoc(K'" \ S;M"); 

(ii) U{x) — )• +00 as X — )• 5. 

(iii) U satisfies 

ess sup ly^^'^jl < 00 VM>0. (36) 

u{x)<M 1 + m 

Theorem 6.1. Under assumptions (i), (ii), (iii), the continuity equation (|22p has existence and 
uniqueness in L~([0,r];L^ n L°°(R'')). 
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Proof. (Uniqueness) Let wt G Lf{[0,T];L^ n L°°(M'^)) be a solution to ([22]), and consider a 
smooth compactly supported function 0:^—7- ]R+. Set E = E{x,p) := + U{x). Then, 

since U is locally Lipschitz on sublevels {U < £} for any £ G M (by (i)-(ii)), (j) o E \s uniformly 
bounded and locally Lipschitz in W^. Moreover 

(V(0 o E) {z),b{z)) = ^'{E{z)){VE{z), b{z)) = for ^'^-a.e. z E 

and we easily deduce that also {(p o E)wt S L'^ {[0,T]; D L°°(M'^)) solves 1^. Let M > be 
large enough so that supp(/> C [— M, M], and let : M — t- M"*" be a smooth cut-off function such 
that V = 1 on [-M, M]. Then (f) o E = o E){(l) o E), which implies that ((/) o E)wt solves (|22|) 
with the vector field b := {ip o E)b. Now, thanks to (i)-(iii), it is easily seen that the following 
properties hold: 

b(^BVxoM'^;^'^), esssup^^^ < oo. (37) 

1 + \z\ 

Indeed, the first one is a direct consequence of (i)-(ii), while the second one follows from (ii)-(iii) 
and the simple estimate 

esssup ' ^ /, < sup ^ , + esssup \ — \^ < oo VM' > 0. 



E{z)<M' 1 + V 1 + IpI / \U{x)<AI' 1 + / 

Thanks to (j37p . we can apply |2l Theorems 34 and 26] to deduce that {(po E)wt is unique, given 
the initial condition ^lq = {(po E)wq^'^. Since E{z) is finite for ^"^-a.e. z, by the arbitrariness 
of (f) we easily obtain that wt is unique, given the initial condition wq. 

(Existence) We now want to prove existence of solutions in L:^([0,T];L^nL°°(]R'^)). Let w € L^n 
L°° iW^) be nonnegative and let us consider a sequence of smooth globally Lipschitz functions 
with I VVfc ~ V^7| —7- in Lj'-^^(M"'); standard results imply the existence of nonnegative solutions 
w'^ to the continuity equation with velocity b^ := (p, — VVfc) with Wq = w, J^dW^dxdp = 
f^aW^dxdp and with ||i(;^||oo < H'W^q IIoo (they are the push forward of Wq under the flow map 
of b^). Since cf) i— )• f^aWf4>dxdp are equi-continuous for all (j) S C^(M'^), we can assume the 
existence of w e L~ ([0, T]; n L°°{W^)) with Wt weakly, in the duality with C^iW^), for 

all t > 0. Taking the limit as A; — )• oo immediately gives that wt is a solution to (|22p . □ 

Theorem 6.2. Under assumptions (i), (ii), (Hi), the u-RLF (^x{t,x,p),p(t,x,p)^ in M^" and 
the u-RLF fj,{t, ^) in ,^(]R^'") relative to b(x,p) := {p, — VJ7(x)) exist and are unique. They are 
related by 

/x(t,/i)=/ \x(t,x,p),p{t,x,p))dl^{x,p)- (38) 

Proof. Existence and uniqueness of the z^-RLF in M*^ follow by Theorem 16.11 and Theorem 14.11 
The uniqueness of the 1/-RLF in I^{W^^ and its relation with the i^-RLF are a consequence 
respectively of Theorem 14.41 and Theorem 14.21 □ 
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7 Estimates on solutions to ([T]) and on error terms 

In this section we collect some a-priori estimates on solutions to ([1]) and on the error terms 
S'e{U,ip)^ S'^{U,ip), appearing respectively in ([9]) and (|12p . 

We recall that the Husimi transform tp i— t- W^V can be defined in terms of convolution of the 
Wigner transform with the 2n-dimensional Gaussian kernel with variance e/2 

-(|x|2+|p|2)/e 

:= = G(")(x)GW(p), (39) 

namely Wsip = (WeV) * ci^"-*- It 

turns out that the asymptotic behaviour as e — t- is the same 
for the Wigner and the Husimi transform (see also (|45p below for a more precise statement). 

For later use, we recall that the x marginal of Weip is the position density Also, the 

change of variables 

h + iy = - ,40) 

[x- = u 

and a simple computation show that the p marginal of Wsip is the momentum density, namely 
(27re)^"| J^^P(p/e)^" (strictly speaking these identities are only true in the sense of principal 
values, since Wsip, despite tending to zero as |(x,p)| — ?• oo, does not in general belong to L^). 
Since the Gaussian kernel G^^\x,p) in (|39p has a product structure, it turns out that 

Ws'il^{x,p)dp= [ \^lJ\'^{x-x')G'^;'^x')dx', (41) 



m\A 



[ WM^,p)dx=(^y [ \Ti;m^^)Gi^\p')dp'. (42) 

Since W^ip is nonnegative (see Section |8] for details) the two identities above hold in the standard 
sense. 

As in [22J we shall consider the completion A of C^(M^") with respect to the norm 

:= / sup \Tp^\{x,y)dy ^ G C,°^(]R2"), (43) 

where J^p denotes the partial Fourier transform with respect to p. It is easily seen that sup \ 
II ^11^, hence A is contained in Cb(M^") and ^(M?^) canonically embeds into A' (the embedding 
is injective by the density of C^°°(M2")). The norm of A is technically convenient because of the 
simple estimate 

/ ^W.^pdxdp <-^y\U\ml (44) 
Since for all if G C^(M^'^) one has if * G^^^^ — ;> (^5 in ^ as e | 0, it follows that 

lim / ipWs^dxdp— / ipWsipdxdp = uniformly on bounded subsets of L^(M'^; C). 

£4-0 JiRd 

(45) 

This will obviously be an ingredient in transferring the dynamical properties from the Wigner 
to the Husimi transforms. 
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7.1 The PDE satisfied by the Husimi transforms 

In this short section we see how (j9]) is modified in passing from the Wigner to the Husimi 
transform. Denoting by T(^y^q) the translation in phase space induced by {y,q) G x M", from 
© we get 

dtr(y,q)Well't + {p-q)- '^xT(y^q)We^Pt = ^y^q)S'e{U,'tpf) 

in the sense of distributions. Since W^ipf is an average of translates of W^ipt: g^t (still in the 
sense of distributions) 

dtWs^Pt + p ■ V^WM = ^t) * Gf""^ + V~eVx ■ [Wei^l * Gi'"^], (46) 

where 

Gf")(y,g):=^Gf")(y,g). (47) 

Indeed, we have 

- / q- Vx^y,q)We^PtG^^''Hy, q) dydq = -^eV^ ■ [W.^t * G^^% 

Although we will not use it here, let us mention that it is possible to derive a closed equation 
(i.e. not involving Wsipt) for Weipt (see |8] and |9], |T^ for applications to the semiclassical limit 
in strong topology). 



7.2 Assumptions on U 

We assume that n = 3M, x = {x\, . . . , xm) G (K^)*'^ and U = Ug + Ub, with Ug the (repulsive) 
Coulomb potential 

l<«<j<A/ ' 

with Zi > 0, and Ub globally bounded and Lipschitz, with VUb G -Bl^oc(K"; 1^")- 
In this context the singular set S of Section [5] and Section [6] is given by 

S = \^ Sij with Sij := {x G M" : Xi = xj} 

l<i<j<M 



Us{x) > (49) 
dist [X, D) 



and therefore 



with c > depending only on the numbers Zi in 

The vector field b = (p, —'VU) satisfies the assumptions (a)-(b) of Section [5] and the assump- 
tions (i)-(iii) of Section El so that the z^-RLF in M^" and the v-RLF in ^(R^") relative to b 
exists and are unique, and the stability result of Section [S] can be applied, as we will show in 
Section (9] 
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7.3 Estimates on solutions to ([T]) 
Conserved quantities. 

A priori estimate. [3 Lemma 5.1]. 

sup/ U]m^dx< [ \Hs^P'o\^dx + 2sup\Ub\{ [{ij'o,Hsi;'o)dx + sup\Ub\). 

teR JR" JR" J 

Tightness in space. jT] Lemma 3.3]. 



sup 

te[-T,T] yR"\B2fl 



(50) 
(51) 

(52) 
(53) 



with c depending only on n. 

7.4 Estimates and convergence of S'^{Ub,'tp) 

In this section we prove estimates and convergence of the term S'i;{Uh,^), as defined in (jlOp . In 
particular we use averaging with respect to the "random" parameter w to derive new estimates 
on ^^{Vjip^), with V Lipschitz only, so that the estimates are applicable to V = Uf,- 

The first basic estimate on S's(V,^p), for tp with unit norm, can be obtained, when V is 
Lipschitz, by estimating the difference quotient in the square brackets in (|10p with the Lipschitz 
constant: 



(oeiV, tp)(l) dx dp 



- 7^r^ll^^ll°o / \y\ \^p4>\{x,y)dy. (54) 

In order to derive a more refined estimate we consider families 'i/'^ indexed by a parameter 
w G W, with {W, P) probability space, satisfying: 



sup sup 



sup sup 

e>0 xeM" JW 



Weijl,{x,p)dF{w) < OO, 

* G\72|^(x) dF{w) < C{X) < OO VA > 0. 



e>0 (a;,p)eR2" JW 
(n) 1 2 



(55) 
(56) 



Under these assumptions, our first convergence result reads as follows: 



Theorem 7.1 (Convergence of error term, I). Let ip^ G L^(]R"';C) be normalized wavefunctions 
satisfying (155 p . (|56p and let V : M" — )■ M 6e Lipschitz. Then 



lim 



W 



^siV,i^'Jcpdxdp+ / {VV,Vp(l))Wsi;l,dxdp 



dF{w) =0 G C^iR^'^). (57) 
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Proof. The proof is achieved by a density argument. The first remark is that linear combinations 
of tensor functions (j){x,p) = (j)i{x)(j)2{p), with 0j G C^(M"), are dense for the norm considered 
in (|54p . In this way, we are led to prove convergence in the case when (j){x,p) = (j)i{x)(j)2{p)- The 
second remark is that convergence surely holds if V is of class (by the arguments in |22) , , 
see also the splitting argument in the y space in the proof of Theorem 17. 3p . Hence, combining 
the two remarks and using the linearity of the error term with respect to the potential V, we can 
prove convergence by a density argument, by approximating V uniformly and in W^''^ topology 
on the support of by potentials £ C^(]R") with uniformly Lipschitz constants; then, setting 
= {V — Vk)(j)i and choosing a sequence in Proposition 17.21 converging slowly to 0, in such 
a way that ||VAfe||2 — )• much faster than l/C(Afc), we obtain 



lim sup 



^siV -Vk,i^'J{x,p)Mx)Mp) dxdp 



dF{w) = 0. 



As for the term in (|57|) involving the Wigner transforms, we can use (|55p to obtain that 



lim sup sup 

fc-5.oo e>OJw 



W,^tn{'^iV-Vk),Vh)h dxdp 
can be estimated from above with a constant multiple of 



dF{w) 



lim sup / \(t>i\\W - Wk\ dx \V(t>2\{p) dp = 0. 

fc-s>oo JR" JR" 



□ 



We shall actually use the conclusion of Theorem 17.11 in the form 



lim 

e->0 



W 



Se{V,ijl,)<P*Gf''Uxdp+ / {VV,Vp<l})Wei^l,dxdp 



u;) = 



V0 G C,°°(M^") 
(58) 

with (f) replaced by G^e^^ in the first summand, in the factor of S^{y,ipl^)\ this formulation is 
equivalent thanks to (|54p . 

Proposition 7.2 (A priori estimate). Let tpf^ G L^(M";C) be unitary wavefunctions satisfying 
(1561) and let (pi, (1)2 ^ Q°°(]R"). Then, for allV ■.M.'^^R Lipschitz and all X>0, we have that 



dF{w) 



(59) 



^eiy,1ptv)ix,p)(l)l{x)(j)2{p) dxdp 

can be estimated from above with 

||0i||oo||Vy||oo / \y\\Tp(P2{y)-J'p(p2*G'^;^\y)\dy + ^/X\\VA\U\Tp(l)2\\i [ \u\G'('\u)du 

+^/C{X)\\VA\\2 [ \z\\Tp(p2\{z)dz+\\V\U\V(l)i\\oo [ \y\\Tp(l)2*G'-;^^\{y)dy. (60) 
Jr" Jr" 

where A := V(pi and C(A) is given in 
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Proof. Set (j)2 = J'p4'2', since (|54p gives that 

can be estimated from above with H^i ||oo ||Vy ||oo / |y||<^2(y) — (i>2 * G^x\y)\ dy we recognize the 
first error term in (|60p and we will estimate the integral of S'eiV, against (j)i{x)(f)2{p)e~^^^'''^^'^ , 
namely 

r I V{x + |y) - V{x - ly) ^^^^^^^ ^ + - '-y)dxdydF{w). 

In addition, we split this expression as the sum of three terms, namely 



II := [ [ Vix+'-y) '^'^''^ ~ '^'^'^ + ^yU 2*G^;:\y)^PU^+'-y)^Ux - '-y) dxdydF{w), (62) 

III :=- [ I V{x- iy)MEl^l±^^, , Gi")(y)C(^ + '-y)^Pf,ix - '::y) dxdydF{w). 

(63) 

The most difficult term to estimate is (j6T|) , since both (|62p and (|63p can be easily estimated from 

above with ^||T^||oo||V0i||oo /ign \y\\4>2 * G^\\{y) dy. We first perform some manipulations of this 
expression, omitting for simplicity the integration w.r.t. w] then we will estimate the resulting 
terms taking (|56p into account. 

We expand the convolution product and make the change of variables (|40p to get 

A{U) - A{u') \ez-{u-u')\ 



Now, the term containing A{u) is equal to 



e i)l,{u)il)l,{u')(l)2{z)dudu' dz. (64) 



- / (ylC) * Gg(u' + ez)ijl,{u')^2{z)du'dz (65) 
and the term containing A{u') is equal to 

- / A(n')C*Gi'!l(n' + ez)V;^.^2(^)(i^i'd^. (66) 
Now, subtract (|66p from (|65p to get that (|64p equals i?^ ^ + 2) where 

^^,1 := - / [(^O * gJ^U^z' + ez) - A{u' + ez)^ * G'J'.U^/' + ez)] V^^(n')02(^)d^^'dz 
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and 



^w,2 '■- 



1 



[A{u' + ez) - * g[%{u' + £z)i:l,{u>)^2{z)du' dz. 



Thus, the apriori estimate on the expression in (|6ip can be achieved by estimating the integrals 
of the error terms ^ w.r.t. w. 
Writing i?^ ^ in the form 



— ^ ^ -G^!^{{u)^l,{u')il:lj{u +ez-u) dudu'dz 

e 



we can estimate from above J"^ \ R^ ^ \ dF{w) by 



IIV^Iloo 

and then by 

VX\\VA 



w 



^G\"i(n)|V'^|(u')|d(u' + ez-u) dudu'dzd¥{w) 

1 vIR" £ 



W JR" 



mz) 



'ne{u)\'iptn\{u')\'ipl,\{u' + ez-u) dudu dzdF{w) 



where ri£{u) := G^"2(^)|^|/(VAe) is a family of convolution kernels uniformly bounded in L} by 

/ \u\G^^\u) du. Using the convolution estimate \\a * %||2 < ||a||2||f?e||i we can finally bound this 

term with \/A||V74||oo||i^2||i / W\G^^\u) du. 

We can estimate from above 2 1 dF{w) using (|56p to get 



\A{u' + ez) - A{u') 



\iljl,\'^{u')dF{w)dudz. 



w 



Then we can use the standard estimate on difference quotients of W^'"^ functions to bound 
this last expression with 



□ 



7C(A)||VA||2 / \zU2\{z)dz. 
This completes the estimate of the term in (|6ip and the proof. 

7.5 Estimates and convergence of S'^(Us,^) 

In the case of the Coulomb potential we follow a specific argument borrowed from [T] proof of 
Theorem l.l(ii)]), based on the inequality 



1 



\z + w/2\ \z-w/2\ 



< 



\w\ 



\z + w/2\\z -w/2\ 



(67) 



with z = {xi — Xj) G M^, w = e{yi — yj) G M^. By estimating the difference quotients of Us as in 
(1671) we obtain: 



S'e{Us,ip)(j)dxdp 



< 



C* / \y\ sup \J'p<j)\{x',y)dy / C/flV'P 



dx, 



(68) 
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with C-t depending only on the numbers Zi in (|48p . 

Now we can state the convergence of S'i;{Us,ip^)] the particular form of the statement, with 
convolution on (f) on one side and convolution on Wsip'^ on the other side (namely the Husimi 
transform), is motivated by the goal we have in mind, namely the fact that the Husimi transforms 
asymptotically satisfy the Liouville equation. 

Theorem 7.3 (Convergence of error term, II). Let ip^ S L^(M"';C) be unitary wavefunctions 
satisfying 

sup / U^l'ip^l'^dx < oo. (69) 

e>0 JR" 

Then 

lim / ^,{Us,^p')(l)*G^J'Uxdp+ [ {VUs,Vp(j))We^lj' dxdp = V0 G C;?°(m2" \ (S x M")). 

(70) 

Proof. First of all, we see that we can apply (|45p with if = (VC/s,Vp(/)) to replace the in- 
tegrals J^2n{^Us,'Vp(f))W£ijj^ dxdp with J^2n{'^Us,Vp(l))Ws'ijj^ dxdp in the verification of (fTOj) . 

Analogously, using (|69p and (|68p we see that we can replace J^2n <^e{Us,ip^)(j) * G^"^ dxdp with 
/]R2n ^eiUs, ip'')4>dxdp. Thus, we are led to show the convergence 

lim / S'siUs,^p'')(pdxdp + / {VUs,Vp^)We^p' dxdp = V(/) G C^{R'^'' \{Sx M")) . (71) 
Since 

S'siUs,'ip^)(l)dxdp = -—^ [ ^'^^ ^i;%x+^y)ij{x - ^y)Tp(l){x,y)dxdy 

[^T^r iR2" e 2 2 

we can split the region of integration in two parts, where y/£\y\ > 1 and where -v/e|y| — 1- The 
contribution of the first region can be estimated as in (|68p . with 

/ \y\ sup \J^pcl)\{x',y)dy / U^\il^'\'^ dx, 

■J{V^\y\>l} x' JR" 



which is infinitesimal, using (|69p again, as e — )• 0. Since 

{VUs{x),y) 



Us{x + ^y)-Us{x-^y) 



uniformly as \/£\y\ < 1 and x belongs to a compact subset of M" \ S*, the contribution of the 
second part is the same as that of 



{VUs{x),y)ip^{x + ^y)V'(2; - ^y)Tp(j){x, y) dxdy 



(27r)'^ J^2n ' """^ ' 2"'^^ 2 
which coincides with 

{VUs,'^p4>)Weil)''{x,p) dxdp. 



□ 
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8 L°°- estimates on averages of t/j 

In this section we consider a family of solutions ipf^ to the Schrodinger equation ([T|) indexed by 
a parameter w, and derive new estimates on their averages. In particular we obtain pointwise 
upper bounds on Husimi transforms. 

One of the main advantages of the Husimi transform is that it is non-negative: indeed, with 
the change of variables (j40)) and simple computations (see |22) for more details), it can be written 
as 

^si^iy^p) = J^{p''^lr>'<Pl,r>) = (^KV','^^,p>P, (72) 
where (-, ■) is the scalar product on L^(M"; C), 

and p't' : L'^{R"-;C) L2(]R";C) is the orthogonal projector onto G L'^{W;C): 



(73) 



[p'^(j)]ix) := / (l)ix')ipix') dx' ^(x). 
Proposition 8.1 (L°° estimates). Let ipf^ £ L'^{M."';C) be satisfying the operator inequalities 

i. f pi'l d¥{w) < Cld Ve > 0. 
e" Jw 

Then: 

(a) for all y G M" and e, X > we have 

C 



w 



(h) for all {y,p) € M^" and e > we have 

Weijl,{y,p)dF{w) < C. 



I 

Jw 



Proof. The proof of (a) follows by applying the uniform operator inequality to the functions 
(26)^^/2 (7rA)"/^G^'^J^2(- - y), whose norm is 1, to get 



e"A"/2 / \^|;f^*G'•^^^,\\y)dF{w)<Ce^. 
Jw 



The proof of (b) is analogous, it is based on (|72p and on the insertion of the functions (f)^^ in 
r3]l in the operator inequality, taking into account that ||0^^p||2 = e""''^. □ 
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The assumption made in Proposition 18.11 is compatible with the famihes of wavefunctions 
given in (fT3]) . i.e. 



C(x)=e-""/Vo(^^)e*("-P°)/^ </>oGC,2(M"), 0<a<l (74) 

with w = {xo,po). Indeed, in this case one can choose W = M^" with the Borel cr-algebra and 
P = p^^", with yO G n see |17j for details. In the extreme case a = 1 no average w.r.t. 
Pq is needed and one can fix it and choose W = M", obtaining convergence for almost all xq, 
so to speak. The other extreme case a = 0, corresponding to concentration in momentum, is 
analogous. 



9 Main convergence result 

In this section we combine the theory developed in Sections |2H6] with the estimates of the Section [7] 
and Section [SI to obtain convergence of the Wigner/Husimi transforms of solutions to (JT]). In 
particular we shall apply Theorem 15.21 

We consider the assumptions on U stated in Section 17.21 and "random" initial data ^ £ 
H'^{W^; C) with unit norm in ([1]) indexed hy w ^ W, where (W, J-, P) is a suitable probability 
space. Denoting by ipf^ the corresponding Schrodinger evolutions, the basic assumptions we 
need for the initial data are 



sup / / \Hi;ipQ ^\'^ dx dF{w) < CO, lim sup 

/ / dx dF{w) = 0; (75) 

e>0 Jw JR" ' -^too e>o Jw Jr"\Br 

— I p^o-^ d¥{w) < CId with C independent of e; (76) 
i{w) := limWeipl^^'^ exists in ^(R'^) for P-a.e. w£W. (77) 

As we discussed in Section |8l (|75p . (|76p . (|77p are compatible with several natural families of 
initial conditions, for instance those described at the end of the introduction, see (|13p or (|74p . 
In addition, the unitary character of the Schrodinger evolution immediately gives 

— / /'V dF{w) < CU Ve > 0, t > 0, (78) 



where C is the same constant as in (|76p . 

We state (|79p below using the Husimi transforms, but (|45p can be used to show convergence 
of Wigner transforms, in the form used in (|16p in the introduction. 

Theorem 9.1. For U as in Section \7. ^ and under assumptions (|75p . (j76p . (|77p . we have 

lim/ sup d^{Wsi^f^,fi{t,iiw)))d¥{w) = 0, (79) 

t6[-T,T] 

for all T >0, where u = i^F e ^(^(M^")) and iJi{t,j2) is the v-RLF in 
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Proof. Our goal is to apply Theorem 15.21 (with a continuous parameter e) and Remark 15.31 with 
ieiw) := Wei^Q u)^'^'^ and fj,^{t,ie{w)) = Weipf^^'^'^. The convergence ([79]) will be a direct 
consequence of (|3ip . We shall work in the time interval [0, T], the proof in the time interval 
[— T, 0] being the same, up to a time reversal. First of all we notice that (|75p and (jSTj) give 

sup sup / / \ H eiif ^\'^ dx (IF (w) < oo. (80) 

£>o teR Jw Jr" 

In particular, by an integration by parts, we have also 

sup sup / / \eV'iljf^fdxdF{w)<oo. (81) 

e>0 teM Jw JR" 

(i) (asymptotic regularity). By (|78p and Proposition I8.ir b) we have the uniform estimate (in e, 
t and {x,p)) 

[ Wsrt^Jx,p)d¥{w)<C. (82) 
Jw 

In particular we have uniform and not only asymptotic regularity, therefore Remark 15.31 applies. 

(ii) (uniform decay away from the singularity). We check (|29p with [5 = 2 and S equal to the 
singular set of C/^, namely 

suplimsup III 27 ^ W si^l^dx dp dt dW {w) < oo. (83) 

5>o e--s>o Jw Jo Jb^ dist {x,S) + 5 

We use (|4ip and the inequality 



dist2(x,5) + 5 ' "dist2(x,5)' 
which holds in Bn for e < e((5, i?) to deduce (|83p from 

limsup / /" / 27 r|V'fu,Pc?2;(ii(iP(w;) < cx). (84) 

£^0 Jw Jo Jk" dist (x, S) 

In turn, this inequality follows by (|52p and (|49p . taking (|75p into account, 
(iii) (space tightness). We have to check that for all 5 > it holds: 

lim p( I It; G 14^ : sup sup / W^eV'fw dxdp> d\] =0. 

V- e>0 te[0,T] yR2"\Bfl ' V 

Considering the cube Cr containing Bji, this tightness property can be checked separately for 
the first and the second marginals of PVeVf using (|4ip . (|42p . it is not hard to see that it suffices 
to check the analogous property for the marginals of the corresponding Wigner transforms; for 
the first marginals, tightness is a direct consequence of (|53p and (|75p . For the second marginals, 
we use and the identity 



\p\'^Wetp dx dp 



1 - ' 



\p\^dp= / \£Vij\^dx (85) 
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with ijj = ipl 

(iv) (time tightness). We need to show that for all (p £ C^(M^") it holds 



lim F{ iweW : 

A/foo 



Wei^l^ dx dp 



dt > M 



uniformly in e. Equivalently, we can consider the limit 



lim ] 

Mtoo 



w eW : 







(t'eWei^lwdxdp 



dt > M 



(86) 



where 



1 * Ge^"^ . According to 



the time derivative in the formula above consists of two 



terms, J{p,'Vx4>e)W£ip^y^dxdp and / ^^iUjipf dx dp and we need only to show a property 
analogous to (Uni) for these two terms. Since (p G C^(]R2"), ||(p,V x</'e)||^ are easily seen to be 
uniformly bounded, hence the first term can be estimated using (|44p . The second term can be 
estimated using (|54p for Ub and (|68p for Ug, taking (|52p and (|75p into account, 
(v) (limit continuity equation). We have to show that 



lim 



cj)We^l^dxdp + ip{t) [ {h.VcfyWei^l^dxdp 



dt 



dF{w) = 



oo(Tn2n 



for all ^ e C, 
validity of the limits 



\{Sx M")), G C^(0,r). Taking (06]) into account, this is implied by the 



lim sup 

^-l-O te[o,T] Jw 



ciP(u;) = 0, (87) 



lim y/e 

£4,0 



w Jo 



<pVx ■ [WMn, * Gi'"^] dxdp 



dtdF{w) = 0. 



Verification of (|87p . We can consider separately the contributions of Ub and Us- For the Ub 
contribution we apply Theorem 17. 1[ in the form stated in (|58p ; the assumptions (|55p and (|56p of 
that theorem are fulfilled in view of (|76p and Proposition 18.11 For the Ug contribution we apply 
(|70p of Theorem 17. 3| the assumption (|69p of that theorem is fulfilled in view of assumption ()75p 
on the initial data and (|52p . ensuring propagation in time. 
Verification of ()88p . This is easy, taking into account the fact that 



{WM^*&^''\V,cP)dxdp 



WM^V,-[cP*G^^''^]dxdp 



are uniformly bounded because G'^^\ defined in ()47p . are uniformly bounded in L^( 



□ 
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